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ff^ Abstract: In this paper the Glauber- Gribov approach for deep-inelastic scattering (DIS) with nuclei is 
T-H developed in N=4 SYM. It is shown that the amplitude displays the same general properties, such as 
^ geometrical scaling, as is the case in the high density QCD approach. We found that the quantum effects 
leading to the graviton reggeization, give rise to an imaginary part of the nucleon amplitude, which makes 
j> jthe DIS in N=4 SYM almost identical to the one expected in high density QCD. We concluded that the 
^ impact parameter dependence of the nucleon amplitude is very essential for N=4 SYM, and the entire 
$H kinematic region can be divided into three regions which are discussed in the paper. We revisited the 
■ - ■ dipole description for DIS and proposed a new renormalized Lagrangian for the shock wave formalism 
which reproduces the Glauber-Gribov approach in a certain kinematic region. However the saturation 
momentum turns out to be independent of energy, as it has been discussed by Albacete, Kovchegov and 
Taliotis. We discuss the physical meaning of such a saturation momentum Qs{A) and argue that one can 
consider only Q > Qs{A) within the shock wave approximation. 
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1. Introduction 

The goal of this paper is very modest and pragmatic: to write a Glauber-type formula for deep inelastic 
scattering (DIS) with a nucleus in N=4 SYM. N=4 SYM at weak couplings is similar to our microscopic 
theory of QCD, with gauge colour group SU(iVc). The high energy amplitude in this theory is given by 
the exchange of the BFKL Pomeron, like in QCD On the other hand, the AdS/CFT correspondence 
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[Q] allows us to calculate this amplitude in the strong coupling limit, where it reveals a Regge behavior 
(see Ref. |3|, ^, |5| and references therein). Therefore, in principle, considering the high energy scattering 
amplitude in N=4 SYM, we can guess what physics phenomena could be important in QCD, in the limit 
of strong coupling. 

The simplest and most informative process that allows to study physics in the region between short 
distances and long distances, is DIS in the wide range of photon virtualities Q. Since the typical distances 
are r oc we can approach the long distance physics at small values of Q. In QCD, we see three different 
regions for DIS: 

1. S> Q1{x) where Ql{x) is the new scale: saturation momentum (see Refs. 0, ^, and a short 

but beautiful review in Ref. @). At such large we can use a linear evolution equation, namely 
the DGLAP equation and the BFKL equation ||l2|, and all advantages of the operator Product 
Expansion |R 



^QCD ^ ^ Qli^)- In this region the density of partons (gluons) is so large that we cannot 
use here the methods of perturbative QCD. However, the QCD couplings are still small here, since 
the typical distances in this kinematic region are r (xl/Qs{x), and Qs{x) ^ ^qcd ■ This fact allows 



us to suggest a theoretical approach in this region, based on non-linear equations [[L4, 15, [L 



< ^QCD- rigorous theoretical approach has been developed in this region in QCD. In high 
energy phenomenology, we describe this region with the soft Pomeron. However, quite a different 
phenomenological approach has been tried in this region, namely, that the saturation scale determines 
the physics inside this domain, and instead of the soft Pomeron, we can use the scattering amplitude 



in the saturation region (see Refs. |17, |l^). Our intention is to use the input from our N=4 SYM 
experience, to penetrate this domain. 

It turns out that N=4 SYM leads to normal QCD like physics in the first region, with OPE and linear 
equations (see Refs. It has been shown in Ref. |[l^ that the DIS densities reach saturation at some 

value of momentum {Qs{x)). However, the physical picture inside the saturation domain turns out to be 
completely different [0] , in the sense that there are no partons in this region and the main contribution 
stems from diffractive processes when the target (proton) either remains intact, or is slightly excited. Such 
a picture not only contradicts the common wisdom, but also contradicts available experimental data. 

In this paper we would like to develop a systematic approach to DIS with a nucleus, based on the 
eikonal formula. In QCD the most reliable approach has been developed for this particular case, since a 
new parameter appears agA^^^ ~ 1, which allows to prove the non-linear equation [^]. 

2. Eikonal approximation for scattering with nuclei. 
2.1 General approach 

It is well known that the eikonal approach is based on two main ideas [pO|, |2l[| . The first one is the fact 
that the value of typical impact parameter for the interaction with a proton is much smaller than the 
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typical impact parameter for the nucleon distribution in a nucleus. Using this idea, we can easily express 
the amphtude for interaction with a nucleus via the interaction amplitude with a nucleon. Indeed, let 
us consider a simple example when the amplitude of interaction with the nucleon is small. Consider for 
example deep inelastic scattering (DIS) with a nucleon. The DIS amplitude for the virtual photon (7*) 
interaction with the nuclear target {A), can be written as follows 

A{-f*A;s,b) = J A{-f*N;s,b') S (b-b'^ J cfb' A {-f* N; s,b') S (b) (2.1) 

where S (6) is the distribution of nucleons in the nucleus, normalized as J d?bS (b) = A, where A is the 
number of nucleons in a nucleus. In Eq. (|2.lD we use the fact that \b — b'\ « Ra ^ Rn ~ . Ra is 
the nucleus radius while Rn is the nucleon size. J cPb' A{'y* N; s,b') is equal to the forward scattering 
amplitude Ai\i{s, t = 0). In the original Glauber-Gribov approach it was assumed that ^Ar(s, t = 0) at high 
energy is mostly imaginary, and ImAN{s,t = 0) = 0"^^* 

The second important observation is the fact that at high energies the longitudinal and transverse 
degrees of freedom are factorized in such a way, that in first approximation the interactions with many 
nucleons in a nucleus will affect the transverse degrees of freedom and the impact parameter distribution, 
but we can neglect the feedback of these interactions on the momentum and the trajectory of the fast 
projectile. In other words, we can use the eikonal approximation for high energy scattering. 

To illustrate this point, let us consider the interaction of the fast particle with the nucleus at rest, as 
it is shown in Fig. ||. 



a) o) 




—V —Pl — p-2 



Figure 1: The single (Fig. |^a ) and double (Fig. |i|-b) rescattering with heavy nucleus. 

First we demonstrate that the momentum transferred q in Fig. |l[ is transverse at high energy. For the 
nuclear target, it is preferable to discuss a process in the rest frame of a nucleus. Describing the nucleus 
in the non-relativistic approach, we consider that the kinetic energy of a nucleon is much smaller than its 
momentum, namely, /2m <C |p| ~ ^/Ra where Ra ^ Rn- Since after rescattering, the nucleon with 

*It should be noticed that such normahzation of the amphtude is a bit unusual for high energy physics since the amplitude, 
calculated from the Feynman diagrams, has a different normalization, namely, ImA — satot- We call the first one as non- 
relativistic while the amplitude of Feynman diagrams will be called relativistic. 
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momentum p — q is still in the same nucleus, qo = Po — {p — q)o = p'^/2m — \p — q\^ /2m <C |g| . In our frame, 
s = 2E Ma where E is the energy of the projectile, and Ma is the mass of the nucleus. At high energy, 
the momentum of the projectile is P = {E, 0, 0, E). Using the fact that = rn^ and {P + g)^ = where 
nip is the mass of projectile, we obtain that 

2P-q = -q^- qo - qz = -q^ /2E (2.2) 
where z is the beam direction. Calculating q^ we have 

= {qo + qz){qo - qz) - ql = -q^/2E{2qo + q^ /2E) -ql (2.3) 

The expression for the diagram of Fig. |l|-a has the following form 



^A(s,g^) = / /f ^4 ■ — — \ — -AN{s,ql;pl,{p-qf); 



— {p — q)^ — ie 



I n w^- r (^^i; tei) 2 \ ■ r (Pi - q; {P^]) (2.4) 

1=1 ' 

where V [pi, {pi}) is the vertex for the transition of the nucleus into A free nucleons. Introducing a new 
variable for the energies of the nucleons, namely, po,i = ^a/A — po^i and noticing that since po,i h^-s the 
interpretation of being the kinetic energy, we anticipate very small values of po,j ^ \Pi\j aiid therefore we 
can neglect Pq ^ . Using this approach, each propagator has the form 

2 2 , 2^ Ma ,^,2 ,^,2 

m -Pi-ie = {--^+^) + '2po,i—^ + \Pi\ - ie = m e + 2po,j "t- + \Pi\ - ie fori < A 

but -p"^ - ie = (- + m^) - 2 ^ po,i + |p|^-ze = me -2 ^ po,i m + - ie (2.5) 

j=i 1=1 

where e = (My^ — Am) /A is the bounding energy per one nucleon in a nucleus, which is much smaller 
than the mass of the lightest hadron. One can see that all propagators for i < A, have poles in po,j in 
the upper semi-plane, while the A-th. propagator has a pole in the lower semi-plane. Closing the contour 
of integration over po,i) on the poles in the lower semi-plane, we obtain the following anticipated result, 
namely 



Aa (s, g^) 



1 



Ae 



2m 



Z^i=2 2m ''^ 



^£ Z^i=l 2m 

r (pi - g; fe}) 



An {s,q\;pl, {p - qf) 



(2.6) 
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The above calculation did not take into account the possible singularities in the nucleon amplitude, 
since their positions are determined by the mass of hadrons po,i ~ ^tt- Closing the contour on these 
singularities, we obtain a smaller contribution of the order of 

Introducing the wave function of the nucleus as follows 



2m 



we can rewrite Eq. (2.5) in the form 



(2.7) 



A 

AA{s,q^;Fig.^-a) = {s,ql) / \\d? ne'^^-''^'^ \^ {{n})\^ {21 

•' i=i 

A 

^ An {s, ql=0) J lid' n e"^^-'''^ I* ({r.}) \^ ^ A^ (s, q\=^) S {ql) 



which is Eq. (2.1) in momentum representation. In deriving Eq. (2.8), we used the fact that in S ((^j^), the 
typical oc 1/Ra-, which is much smaller than the characteristic q±^ in the nucleon amplitude, and which 
can be considered to be a constant as far as the q_\_ dependence is concerned. Now we want to show that 
the diagram of Fig. |^b leads to the following contribution 



AA{s,h-Fig.^-b) = 



cfib'AN{s,b')] S^b) 



(2.9) 



It turns out that Eq. (|2.9| ) can be obtained with the additional assumption that the wave function can be 
factorized as 



^{{ri}) = Yl "^{ri) , which gives 5(6)= / dz|*(6,z)|2 , with f = (6j_, z). (2.10) 
1=1 

This means that there are no correlations between different nucleons in a nucleus. In other words, we 
describe the nucleus as the nucleons that are moving in the external potential in the spirit of the Hartree- 
Fock approach. 

The amplitude for the diagram of Fig. Q-b has the form 
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A^[s^q'■F^g.^-h) = j Jp ^ j H ^ fe}) 



AN{s,k\;pl,{pi-kf)] 



w? — P2 — ie frfi — pfl — ie 



An {s,{q-k)'{;pl,{p2-q + kf) 



m? — {pi — k)"^ — ie 
X 7 — TTT2 2 — ~2 -'^iPi-k,P2-q + k,{pi}) (2.11) 

We integrate first over the momentum k. Rewriting d'^k as dkod^ko — kz)d'^k, and closing the contour 
of integration over the variable ko — kz, on the pole {P + fc)^ = rrip, leads to a factor of 27ri/Po. For the 
integration over /cq, we can also close the contour on one of the poles: {pi — k)"^ = m? or {p2 — q + k)"^ = m? , 
which can be rewritten as me + 2m (po,i — ^0) ~ (Pi — k)'^ — ie = and me + 2m (po,2 — Qo + ^o) ~ iP2 — Q + k)'^. 
This integration brings an additional factor of 27ri/2m. Therefore, the integration over k leads to the 
following contribution, namely i(f'k/{{2-K)'^ s). Evaluating all the integrations over po,i in the same way as 
we did when calculating the diagram of Fig. |l]-a, we reduce Eq. ( |2.1l| ) to the following expression 



A^(.,,i,F^a.%-h) = '- n |y r (pi ; {p.}) ^ ^ _ ^ ' _ A„ (s,kl:pl) (2.12) 

X A, (, - k)l;pi) r (p. - {«}) 

2m 2m Z^i=3 2m 



Eq. ( |2.12 ) can be easily rewritten in coordinate representation, by introducing the wave function of 
Eq. (|2.7D , namely 



Aa {s, h- Fig. I - 6) = -^Al (s, = j dzi J ' dz2 \^ {, 6, zi;b, Z2; {rj) \^ (2.13) 

Using the non-relativistic normalization for the scattering amplitude ( Anr = A/s)^ and Eq. ( p.lOD , 
we can see that we obtain Eq. (|2.g|). It should be noted that the factor 1/2 stems from the Z2 integration, 
which is not restricted in Eq. (^.9|), in contrast with Eq. ( |2.13D . All calculations above have been done to 
illustrate two points, namely that we do not need to assume that the nucleon amplitude should be pure 
imaginary, but we need to assume a very simple model for the nuclei. 

Calculating the amplitude for the interaction with any number of nucleons in a nucleus, we obtain the 
simple formula for the nucleus scattering amplitude (see a more detailed derivation in Ref. p^), namely. 



^Starting from this equation we will use the notation An and Aa for the non-relativistically normalized amplitudes, hoping 
that it will not lead to any misunderstanding. 
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AA{s,b) = - exp(^i j db' AN{s,b') S {b)j^ (2.14) 

In deriving Eq. ( ^Tll ) we considered the propagators of the projectile and the target (nucleons in a nucleus) 
in flat space but not in AdS^. In the next section we will comment on this but the main argument is very 
simple: the trajectory of a fast moving particle can be replaced by the straight line in curved space as well 
as in flat one. The second assumption was that we considered in Fig. |^b the projectile as the intermediate 
state. 

^ Using the AdS/CFT correspondence we can esti- 

mate the accuracy of this (eikonal) approach in the 
^ N=4 SYM case. Indeed, at first sight we can expect 

from the AdS/CFT correspondence, that the main con- 
^ tribution will stem from the fan diagrams, the first of 

a,) b) which is shown in Fig. ^a, as it happens in this theory 

m in the region of small coupling constant. In fact, from 
Figure 2: The first fan diagram (Fig. &a) for the 

/ -J \ 1 the region of small coupling we expect that (i) this di- 

mteraction of Pomerons (reggeized gravitonsj and the r- o r- \ oa 

eikonal diagram (Fig. |-b). agiam has the contribution of the order of (af /A) s , 

where A oc is the intercept of the BFKL Pomeron; 
(ii) the typical value of y — Yi ~ 1 /A ^ 1 and (iii) the value of this contribution is closely related to 
the process of diffractive dissociation of the projectile. Since Y — Yi ^ 1 it is reasonable to consider the 
exchange of the BFKL Pomeron. The eikonal diagram of Fig. ^b has the same order of magnitude but it 
turns out (see Ref . [p4| ) that this diagram is included in the diagram of Fig. |2|-a in the region of integration 
Y — Yi ~ 1 where we cannot use Pomeron exchange. Therefore, in the weak coupling limit the full set 
of diagrams at high energy can be reduced to the "fan" diagrams. It is worth mentioning that in the 
weak coupling limit the eikonal diagram of Fig. l-b has the same intermediate state as the initial one (the 
colourless dipole) since it turns out that colourless dipoles are diagonalized by the interaction matrix (see 
Ref. PI). 



In the strong coupling limit of N=4 SYM, due to the AdS/CFT correspondence, the strong interaction 
of Pomerons is replaced by the weak interaction of the reggeized gravitons, with intercepts A = 1 — 2/-v/A, 
and therefore in the triple Pomeron diagram the typical value of F - Yi « 1/(A = 1 - 2/^/X) « 1. It 
means that diffraction production, which can contribute and was neglected in the eikonal (Glauber-Gribov) 
approach , is the process in which low masses are produced. For Y — Yi « 1 there are no reasons to replace 
the amplitude by the reggeized graviton exchange. Using the AdS/CFT correspondence we expect that 
in the diagram of Fig. l-b the same as the initial state is produced. On the other hand, the process of 
diffraction production of low mass can be easily taken into account in the eikonal approach, and does 
not change neither the energy nor the impact parameter dependence that has been discussed here. The 
cross section of the diffraction dissociation is proportional to the imaginary part of the reggeized graviton 
exchange which is small of the order of Therefore, at least within this accuracy ( 2/^/X), the 

exchange of two gravitons between the projectile and the target (eikonal diagram of Fig. ^b) prevails. 
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2.2 Nucleon amplitude in N=4 SYM 



The main contribution to the scattering amplitude at high energy in N=4 SYM, stems from the exchange 
of the graviton-l- . The formula for this exchange has been written in Ref.|Q, (see also Ref. [10| for its 
interpretation). In flat space this amplitude has the following form 

Ac,{s,q) (X T^u{pi,P2)G^,yf,'^i (q) T^'jyi {pi,p2) ^—^s^— (2.15) 

where T^^jj is the energy-momentum tensor, and G is the propagator of the massless graviton. The 
last expression in Eq. ( ^.15 ), stems from the fact that for high energies, T^^y = pi^^pi^u, and = —q\ (see 



the previous section). However, we are interested in N=4 SYM in a space with curvature, namely AdS^. 
AdSd+i corresponds to an hyperboloid in d + 2 flat space, namely 



-Y.\ + + Y1 = (2.16) 

7=1 

with curvature R = —d{d — 1)L^. Introducing new coordinates 



Yo + Y-i Yo + Y. 

we reduce the introduced metric to the following form 



(2.17) 



In the flat d + 2 dimensional space, the scalar propagator is the following ( with y+ = Yq + ^-i and 
Y_ = Yo-Y_i) 



G(Y„Y+,Y_;Y,',Y' Y') = /" TT — ^ 

^ ii 2^ (2^)2 ~ *P+P- -ik-X - i-p+Y- - t-p-Y+ 

= (2^)-'^/2-i r dtt-"/^-^ e~^/' (27r)-'^/2-i d^(e)'^/2"i e"?" 



= {2TT)"'^/'^~'^T{d/2) u~'2<^ (2.19) 
In Eq. ( |2.19| ), n is a new variable which is equal to 



* Actually, the graviton in this theory is reggeized but it is easy to take this effect into account (see Refs. js], |5|, p^) 
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„ = + (2,20) 

2z z' 

In Eq. ( |2.19D , we re-write the integration measure of the momenta in d + 2 dimensional space, namely 
nf=i '^Vii ^-s nf=i '^^i dp^dp-, where p-^- and p_ are the conjugated momenta to y_ and y+, respectively. 

However, the propagator of Eq. ( p. 191 ) does not satisfy the correct boundary condition, for example, 
G (X., y+, y-; X;, n, r ) should approach 5 [x - X') as z ^ which is not the case for this equation. 

One of the reasons why this happens, is that we have to guarantee that 1+ > 0^. The easiest way to impose 
such a condition, is to change Eq. ( |2.19| ) to 



G {Xi, Yj^,Y^;X[, y|_, yL 



ndki dp^dp^ 
9^ 



L2vr (2vr)2 Y^U^^ + P+P- P~ 



fyrdh dp+ J_ 

J 27r (27r)2 k"^ 



i ^ Y+ 



~ik-X-i\p+Y- 



(2.21) 



One can see from Eq. (Q) that (J^^q d'^/d'^Xi + d'^/d'^Yo - d'^/d'^Y^i) G {Xi, Y+,Y_;Xl, y|, Y^) is de- 
fined only for 1+ > 0. Therefore, the solution of the equation for the Green's function also will be 
determined only for Y^ > 0. 

Notice that the mass of the graviton is equal to zero even in the AdSd+i space with curvature. Having 
this in mind, the easiest way to find the correct propagator, is to write the wave equation directly in the 
AdSd+i space, assuming that the mass of the graviton is equal to zero, and that G [Xi,Y^,Y_; X'-, Y^,Y!_^ 
is a function of the variable u of Eq. ( 2.2C| ). The action for such a particle has the following form 



5[0] 



(2.22) 



where the metric is given by Eq. ( |2.18| ) . Using Eq. ( |2.22D and Eq. ( ^.18| ), it is easy to obtain the wave 
equation for G {Xi,Y+,Y^; Xl,Y[,YL) = G{u). It has the form ||, || 



-^d^^g>'-'d,G{u) = 0; 



z^VlG{u) + z^+' 



d_ 

dz 



,1 dG{u) 
dz 



0: 



u (n + 2) (n) + {d + l)Guiu) = 0; 



(2.23) 

(2.24) 
(2.25) 



The solution to Eq. ( 2.25 ), which satisfies all the necessary boundary conditions: G (u) 



and 



G(u) 



6{x- x') has the form|22|, |3|, § 



^We thank Chung-I Tan for the fruitful discussion of all aspects of high energy scattering in N=4 SYM, in particular, the 
Y+ > condition. 
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As has been discussed (see Eq. ( |2.15| )), we need an expression for the propagator of the graviton, which 
at high energy depends only on the transverse coordinates for the scattering. Therefore, we need G (u) for 
AdS2+i, which is equal to 

Gsiu) = ^- (2.27) 

with 



2zz' 

where b is the impact parameter for the scattering amplitude. 

For the eikonal formula, we need to evaluate the integral over b, which can be easily done noticing that 

f , . 1 \l + u + + 

d{l + n + Vu{u + 2)\/db^ = -—j^ . , ^, ^ (2.29) 



The result is 



G{z,z') 



I 



(2.30) 



d'bG,{u) = ^ 



Z^ 



{l + u(6 = 0) + ^u{b = Q){u{b = Q) + 2)y 



Z z 



{z^ + z'^ + \z^ — Z' 



I2\\2 



This equation provides us with the factor which enters into Eq. ( p. IS] ), instead of It turns out 

that in curved space we need to change H 



zz' s 



(2.31) 



For calculating the nucleon amplitude, we need to multiply Eq. (2.30) by the coupling constant, and 
integrate over the nucleon wave function Therefore, the nucleon amplitude is equal to 



/ 



dHANis,b) 



J dz' zz' Giz,z')mz')\'^ = igls j dz' \<^> {z') 

dz'|$(z')|2 



z2 z'2 



> 7— Z 



(z2 + z'2 + |z2 - z'2|)2 

iglNcS 



(2.32) 
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Here Qq is the dimensionless constant, which is equal to k|/2L^, where K5 is the five dimensional gravity. 
5q oc where Nc is the number of colours. We do not know anything about the nucleon wave 

function, except that the integral over z' converges, and it is proportional to Nc- Therefore, the amplitude 
is proportional to A]\f oc s/Nc and it is small for sz'^ < ■ It grows and becomes of the order of 1 
due to the reggeization of the graviton. The graviton propagator in Eq. (|2.32| ) should be replaced by the 
propagator of the Pomeron, in the way as has been suggested in Refs. |^ |l^, ^. This modification for our 
case is described in section 5. In Eq. ( 2.32 ), we consider J dz'\^z')'^ = Nc ■ 

As has been discussed, we use the propagator for a fast moving particle in the form 

1 



G [k+,k_;bi - 62; 2^1 



Z2 



5(2) (^h^ - 62) 6 (zi - Z2) 



(2.33) 



Eq. ( |2.33| ) follows directly from Eq. (2.21) . Indeed for large the pole in the integrant of Eq. ( p. 211 ) is 
located at = {k']^—p^p^—ie)/k^ —>■ — ie. Therefore, ^^^-^ A;? can be replaced by /c+(A;_ +ze). 
Substituting this expression in Eq. ( |2.2lD one can see that G (jt^,k^;bi — 62; 2:1 — Z2^ has the form of 
Eq. ( |2.33| ) with an additional factor @{zi + Z2) which is equal to 1. 

Eq. (|3|) for G(k+,k^]bi-h2]zi- z-^ can be derived directly from Eq. (|]2|) and Eq. { ^^ . 
Indeed, going to Fourier transform for coordinates Xi{i = 1, . . . ,d) and to Laplace transform for coordinate 
z we can rewrite Eq. ( 2.24| ) in the form 



k''G'^{{h]-p) - {d-l)pG{h}-p) 
The solution to this equation has the form 



p'G{{h}-p) 



(2.34) 



G{{h]:p) 



k' 



For large A:+ Eq. ( p. 351) leads to 



d-l 
2 



Gi{h};p) 



k+ k^ 



k+ <c {fcx and p} 



k^, 



kik^ 



kj^k- 



k', 



d-l 
2 



1 



fc+ (A;_ — if) 



(2.35) 



(2.36) 



Eq. ( |2.36| ) gives Eq. ( |2.33| ) which we use in our calculations. 



2.3 Eikonal formula in N=4 SYM 

Eq. ( |2.14 ) can be easily rewritten for the case of N=4 SYM in the following way using Eq. ( 2.32 ) 



dz \^p{z) 



6 4 



2^ S{h) 



(2.37) 



where <I>p is the wave function of the projectile. This formula is almost the same as the eikonal formula 
for the hadron-nucleus interaction, except that the nucleon amplitude is purely real in our case. 
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The scattering amplitude at fixed z 



Aa {s,b;z) = i < 1 



can be rewritten in the following way: 



Aa (s, b; z) 



sm 



+ i 2sin^ 



5 (6) 



(2.38) 



(2.39) 



One can see that the real and imaginary part of the amplitude are of the same order in contrast with 
the black disc behavior, for which only the imaginary part survives at high energy. One can see that the 
amplitude of Eq. ( |2.39| ) satisfies the following unitarity constraint 



2ImAA{s,b;z) = \Aa {s,b; z) 
Comparing Eq. ( p. 40 ) with the general unitarity constraint, namely, 

2ImA{s,b;z) = \A {s,b; z) {"^ + Ginei {s,b; z) 



(2.40) 



one can see that Eq. ( p. 39 ) leads to only elastic scattering at high energy, in direct contradiction with our 
intuition based on the parton approach. 

For the general formula of Eq. ( p. 37 ), Eq. ( |2.40| ) means that 



cjtoi = 2 d^b / dz\^p{z)\^ Re< 1 - e'' — 



90 ,4 



Z^S{b) 



o-diff + (^ei = d b / dz\^p{z 



1 _ e*^^ ^'^W 



(2.41) 



In other words, only the processes of diffractive dissociation contribute at high energy. 



3. DIS with nuclei: general formulae 



For calculating DIS, we need to specify the wave function of the projectile in Eq. ([2^. In N=4 SYM 



the natural probe for DIS is 7^-current (7^-boson) |19], and the wave function for this probe satisfies 
Eq. ( |2.24| ). However, in DIS we fix the virtuality of the probe (see Fig. P). It means that in terms of 
Eq. (|I|), ^^ti kf = Therefore, the wave function is described by Eq. (|2.24D with d = 0, but with 



-Q^, and the equation can be rewritten in the form |19] 



-z^Q^^n {Q^z)+z 



d^n{Q\z) ^d^^niQ^z) 



dz 



+ z 



{dz) 



(3.1) 
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Figure 3: DIS with the nuclear target. The wave hne denotes the TZ current (7?.-boson), while the zigzag lines show 
the graviton exchanges. is the virtuality of the probe. 



The solution to Eq. (|3.1|) is 



"^TliQ^z) = QzKiiQz) 



(3.2) 



However, TZ - boson is a vector with d + 1 components. The careful analysis of ref. |1£] shows 
that Eq. ( |3.2[ ) describes only d components of this vector, while the {d + l)-th component has a different 
dependence on Qz. Finally|10|, 



{KfiQz) + KliQz)) 



(3.3) 



The deep inelastic structure function has the following form[^, 



F2{Q\x = Qys) = (3.4) 



Ca'Q^ j d^b J dzz^ {Kf {Q z) + {Q z)) 2i?e|l - exp (^i 



4 X 



z" S (b) 



where C is a dimensionless constant. 

Changing the variable z to C = Q z, one can see that F2 can be written in the form 

F2{Q^x = Q^/s) = C J d%<^{T{Q,x,b)) = (3.5) 
= CQ^ J dHj dCe{KfiO + K^.iO) i?e|l - exp(^iic')| 



where 



_ Q^x _ Q2 

One can see that the DIS structure function shows the geometrical scahng behavior with the saturation 
momentum, which we can find from the equation with r = 1 . It is equal to 



Ql{x) = glN,S{b)/{Ax) « (3.7) 

i Vc X 

Therefore, F2 shows the same main features as F2 in high density QCD 0, ^, namely the 

geometrical scaling behavior, large values of the saturation scale in the region of low x, and the expected 
dependence of A^/'^. Actually, our analysis of Qs repeats the one in Ref. |1C], and the difference 

between them stems from our integration over the impact parameters. 

One can see from Fig. ^ that the function $ has the same behavior as we expected from high density 
QCD, namely it approaches unity at small values of r. Such a behavior looks strange, especially if we 



compare this function with Eq. (|2.39D , which leads to an amplitude that oscillates between and 2. Let 
us consider r > 1. In this case, we can replace the modified Bessel functions (McDonald functions) in 
Eq. (|3.5| ) by their asymptotic expression, namely, Kn{C,) — > y27r/C 6xp(— C), and in this case Eq. (|3.5|) has 
the form 

F2{Q^,x = Q^/s) = C J cfb^iT{Q,x,b)) = (3.8) 
= C2ttQ^ J (fb J dCC^e-^ iJejl - exp (^i i | 

The second term in {. . . } can be estimated by the saddle point method. One can see that the saddle 
point value for = (^sp = (~^''"/3)^''^, and the integral has the following form 



$,p(r) = 1 - y^4pe-(^/^)(-/^)^^^ -> 1 (3.9) 

One can see that at r ^ 0, the exponent e"'^^/^-*^*'^/^^^'^^ 1, but the pre-exponential factor tx r^/^ 
vanishes. However, since (sp !> small values of r we have to use the expression for the modified 
Bessel function at C ^ 0, namely i^n(C) ^ — ~^ 1/C"- Doing this, one can see that C ~ ^"^^^ contributes to 
the integral leading to the behavior of the second term in Eq. ( |3.9D proportional to ^/t. 

The above discussion shows that predictions of high density QCD differ from those of N=4 SYM, only 
in the way that $(r) approaches unity, namely $ — 1 oc exp (-Cln2(l/r)) in high density QCD, and 
•3> — 1 cx exp (— iln(l/r)), in our approach. 

We need to integrate <I> (t(6)) over b (see Eq. ( |3.5|) ), to obtain the total cross section for DIS 
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atot{DIS) = —F2{Q^x = Qys) (3.10) 



lr{b=0) T {b{T)) 

^ C^R\ r ^ ^(-) -here R{r) = #^ 

where r = Tmax is the position of the maximum of the function $(t). The exphcit form of the function 
i?(r) depends on the dependence of S{b) on the impact parameter. We hst below this function for several 
nucleus models: 



R{t) 



T cylindrical nucleus S{b) = {A/ttR\) © {Ra — b) ; 

1 Gaussian form 5(6) = {A/'kR\) exp [-b'^/R\) ; (3-11) 



t(6 = 0)/t^ spherical drop nucleus S (6) = {SA/Att'^ RI 



A) 



62; 



0.8 



0.6 



0.4 



0.2 




Unfortunately, in a realistic model of the nu- 
cleus with the Wood-Saxon form for the 6 depen- 
dence, we cannot give a simple analytical form of 
the function R{t). In Fig. |5|we plot the integral 
over r in Eq. ( p. 101 ), for Gaussian 6 distribution. 
This distribution, being oversimplified, leads to 
correct estimates for the average characteristics 
of nuclei. 



^From Eq. ( 3.10| ), one can see that the total 
cross section for DIS will be 2-kR?^ x ln(l/r(6 = 
0)), once more in accordance with our expecta- 
tion from high density QCD for such 5(6). In 

the case of the Wood-Saxon parameterization, 

b>RA 



2 

ln(l/x) 



5(6) 



exp(— 6//i) which leads to a tot oc 
0). This behavior coincides with the 



Figure 4: The r dependence of function 



expectation of high density QCD. 

Therefore, the Glauber-Gribov approach leads 
to a behavior of the DIS structure function, which fully supports the high density QCD picture, repro- 
ducing the geometrical scaling behavior, and the existence of only one new scale, namely the saturation 
momentum. 
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The main difference between N=4 SYM and high density QCD, hes only in the relation between the 
total cross section and the cross section of diffractive dissociation. That is, atot (DIS) = (Jdiff (DIS) for 
N=4 SYM, and atot {DIS) / adiff {DIS) but 

fdi// > 2^*ot for high density QCD. In N=4 SYM, this equality means that the elastic cross section is 

equal to zero, in sharp contradiction with QCD and any parton interpretation of high energy scattering. 
However, this is a direct consequence of the fact that the graviton has spin 2. Actually, it has been shown 
in Ref.[^ that its spin in N-4 SYM is not exactly 2, but rather jgraviton = jo = 2 — Because of this, 

the amplitude of the interaction with the nucleon is not purely real, as it is given by Eq. but it has 

an imaginary part which is proportional to 2 — jo . Fig. ^ illustrates how this imaginary part influences the 
total and inelastic cross sections. We introduce the functions ^tot and as 

at. = jdH^todr) and a..^ = j dH 

The functions ^tot and are shown in Fig. ^, for the imaginary part of the graviton exchange, which 
is 10% of the real part of the amplitude. One can see that such a small imaginary part generates a large 
inelastic cross section, and therefore the DIS structure function in N=4 SYM, with reggeized graviton, 
leads to a qualitative picture which is very difficult to differentiate from the high density QCD predictions. 




2 4 

hi(l/x(b=0)) 



Figure 5: The integral over r in Eq. ( 3.1C ) as a 
function of T{h = 0) for Gaussian dependence of the 
nucleon density in nucleus versus the impact param- 
eter. 
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Figure 6: The behavior of the total and inelastic 
cross sections for the graviton exchange with 10% 
imaginary part of the amplitude. 



To complete the proof of Eq. ( |3.8D , we need to discuss the contributions from multi-graviton exchanges 
in the nucleon amplitude. At first sight, they should be essential, since each graviton exchange brings in a 
factor (see Eq. 
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A%{s,b) = igos J dz' mz')\'^ z z' Gsiu) Sig^s z^ j dz' \^{z')\'^ z'^ / (b^f (3.12) 



4N1/3 



^From Eq. ( p^ , one can see that the aniphtude j4^(s,6) ^ 1 for b'^ — ^ (^'^ 

z^y\ This means 



that we need to take into account all terms of the order of (^^)". Using the eikonal formula for summing 
such terms, we see that for the nucleon amplitude we have the following expression, instead of the simple 



formula of Eq. 



d^bAN{s,h) = (3.13) 
j dH{l - exp(A^(s,6))} = j d^ftjl - exp^iflfos j dz' \^{z')\'^ z z' G^{u) 



1 In 

The integral over b can be estimated as / d'^bA^ (s,6) oc C tt^q cx (isz ^) ' . Using Eq. (|3l^ ), we 
can rewrite Eq. ( 2.37| ) in the form 



AA{s,b) = i J dz|$p(z)|2|l - e'^'^^'oSib)^ (3.14) 



For DIS we have 



i^2((9^x = QVs) = (3.15) 
Ca'Q^ j d% j dzz^ [kI{Qz) + kI{Qz)) 2Re{l - exp(i7rC" blz*S{b))] 

where C and C are dimensionless constants, whose values are irrelevant for our discussion. 

Performing the integral over z using the asymptotic behavior of modified Bessel functions and the 
saddle point approach, one can see that the saddle point value of z = z^p is equal to 

= iSm " ^ ^'-''^ 

where A is the number of the nucleons in a nucleus. The value of F2 in the saddle point is 

F2 oc exp (-i const Q^/(sA)) (3.17) 

This formula, if correct, leads to a saturation scale oc A/x, in drastic contradiction with the 
prediction of high density QCD, both in the A and s dependencies. However, if we come back to Eq. ( pl|) , 
we obtain 
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A%{s,b) = Sig^szi^P j dz'\^{z')\^ z'^l [lo^f cx ^^(^) < 1 foi" Q^(x ^ 0) (3.18) 

^From these estimates we conclude that the multi-graviton exchange does not contribute to DIS with 
a nuclear target, at low x. 

4. DIS with nuclei: ultra high energy limit 

The result of the previous section is, however, valid only for a limited range of energy. Indeed, we observe 



[IS z 



that the value of the typical impact parameters in the nucleon scattering amplitude ( 6^ = 6q 
grows with energy, and for energies larger than the energy (s = Scrit) when 6o > Ra^ we cannot use the 
eikonal formula in the form of Eq. ( |2.14| ). Indeed, for such large energies, the main assumption of the 
Glauber- Gribov approach does not work. This assumption has been discussed in Eq. ( |2.8| ), which can be 
rewritten in the following way in the case of one graviton exchange 

AA{s,b) = J (fb' An {s,b') S (b-b'^ j dH' An {s,b') S (b) (4.1) 

In Eq. {iA), we assume that in the interaction with one nucleon, the typical impact parameters are 
much smaller than Ra, which gives the scale for the impact parameter distribution in the nuclei. If the 
typical b in the nucleon interaction is larger than Ra, we have to use a different approximation, namely we 
need to rewrite Eq. ( [4.1| ) in the form 

AAis,b) = j (fb' An {s,b') S (b-b'^ ANis,b) j dH' S (b) = A An {s,b) (4.2) 



This equation leads to a new formula for the scattering amplitude with a nucleus, instead of Eq. ( p. 14 ), 
namely, 

AAis,b) = i{l - exp{iAANis,b))) (4.3) 
which leads to an expression for the DIS structure function in the form 

F2{Q\x = Qys) = (4.4) 
CaQ^ j dH j dzz^ {KfiQz) + (Q z)) 2 Re{l - exp {i A An (s , b))} 



where An {s,b) is given by Eq. ( 3.12 ). An can be rewritten at large b, using Eq. (|2.3D and Eq. ( |3.12|) , in 
the form 



- 18 - 



Am {s,b; z) 



I sm 



9^0 Nc 



2\3 



+ 2sin^ 



16 



2\3 



(4.5) 



Substituting Eq. (4^), we do the integral over z using the steepest decent method. The most important 
part of the nucleon amphtude is the imaginary part, which leads to a damping of the interaction matrix 
(5-matrix) at high energies, provided the amplitude tends to unity. The saddle point for z is equal to 



ZSP 



4As cos 



Q 



zV(62)S 



(4.6) 



Taking the integral using the steepest decent method we obtain 



F2 {Q\x = Q'/s) ^ I d'bzlp J^^^^ exp ("-5/4 '"j (4.7) 



where we replaced sines and cosines in Eq. (^]^) and Eq. (^]^) , by unity since these functions cannot change 
the energy and Q dependence of the resulting amplitude. 

^From Eq. (|4.7D , one can see that the typical values of the impact parameters are large and equal to 



hi = ^/{hQzsp) = ^ (-^)' » 4p (4.8) 



The resulting answer for F2 is the following 



F2 = oc a'Q^ = a'Q^Ahx-h (4.9) 

Therefore, we see that we expect a very strange behavior from the point of view of high density QCD, 
both as function of A and x. The origin is clear. N=4 SYM has a massless particle, namely the graviton, 
and because of this the nucleon amplitude falls at large 6^ + z'^, as a power of 1/(6^)^ . Such a 

power-like behavior leads to a typical h which grows as a power of energy, (see Ref.|Q for details), as has 
been demonstrated above. However, as has been shown in Refs. [^, |5[, actually the graviton has a mass 
which is not equal to zero if we dealing with the propagation of the graviton in AdS^. This mass leads to a 
reggeization of the graviton, which has spin jo = 2 — 2/\/A < 2, in the scattering kinematic region where 
the square of the momentum transferred t is negative {t < 0). The fact that there is no massless particle 
in the curved space means that at large b, the amplitude should falls exponentially leading to a log energy 
dependence of the cross section. This is the reason why in the next section we will discuss the exchange of 
the reggeized graviton, and the Glauber- type formula which such an interaction induces. 
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5. DIS with nuclei: graviton reggeization. 



As has been discussed in Ref. for the exchange of the reggeized graviton we need to replace Eq. ( p. 32 ) 
by a more general expression, namely 

/ ,HA.(s.,, . {-/ |1 (ii-tJ^) / ,HG,iu„)] (5.1) 

where 



I (l + n+ y^m('U + 2) 



2-A+{i) 



(5.2) 



with 



A+OO = 2 + a/4 + 2VA(j - 2) = 2 + J2^{j 



Jo) 



Using the definition of u given in Eq. ( ^^281 ) and Eq. (|2]29|), we can easily evaluate the integral over b 
in Eq. (^T^) with the following result 



(fb G3 {u,j) = zz 



2 - A+ (j) 



(1 + u{b = 0) + ^u{b = ^) (n(6 = 0) + 2)) 



2-A+(j) 



(5.3) 



The integral over j in Eq. ( |5.1| ) is a contour integral, and the contour is located to the right of all 
singularities of / G3 (n, j), but to the left of j = 2, and the contour can be drawn to be parallel to 
the imaginary axis. In Eq. ( |5.3D , one can see that our singularity in j stems from the zero of the factor 

z/, we can rewrite the contribution of the square root singularity at 



2 - A+ (j). Denoting -^2^ |j - jo 
j = jo in the following way 



/ 



d^bAt^{s,h) = gl2zz' [cot^ + i] s^°-^ 



(5.4) 



z>z' 



ie+00 



dv 



vr 



exp (^-i/V(2\^) + ii/ In 1 1 + n(6 = 0) + Vu(6 = 0) (u(6 = 0) + 2)|) 



0^2 z z ( cot + I \ s- 



dv 



exp ( —1/ /(2\/A) + «zy In ( — 



In the course of deriving Eq. ( ^.4[) , we neglected in the signature factor the contribution of the term 
j(2\f\)^ considering it to be smaller than jo (jo > v'^|{2^J\). The integral in Eq. (U) can be evaluated 
such that it reduces to the following expression 
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bAj^f {s, b) = g^lz z [ cot + i \ s- 



TT ^/\ In S 



exp 



21ns 



(5.5) 



The result of Eq. ( |5.5D is obtained assuming that A is fixed, but s ^ oo. From Eq. ( p.lD , Eq. (|5^ 
and Eq. (|5.3D , we can recover a diff'erent limit, namely A ^ oo when s ^ 1. Indeed, in this limit 

— > 4 + (\/A/2) (j — 2). Since 2 — A+(j) (\/A/2) (j — 2), we can close the contour on the pole which 
stems from 2 — A+(j) = 0. The signature factor can be rewritten in the form 



cot + I 



Collecting everything together we obtain 



A^oo^ \/A 



(5.6) 



j dHAN{s,b) = gl2zz'\l - i-^l j dHGs{u,2) = gl2zz' \l - j dHG^{u) 

(5.7) 

We have used Eq. (|5.7| ) in our estimates of the value of the imaginary part of the nucleon scattering 
amplitude. Eq. ( ^.71) leads to the exchange of the graviton with a small imaginary part, and this case has 
been considered in detail in this paper. 



We concentrate our efforts on the limit s — > oo,A = Const. For this region we need to use Eq. ( |5.4D for 
the nucleon amplitude. However, even more important for the high energy behavior of the amplitude, is 
the fact that the graviton has a mass in curved space (see Fig. ^). Therefore, the graviton trajectory which 
gives the dependence of the spin of the graviton on its mass, has the intercept jo = o-graviton (0) = 2-2/y/X, 
and the mass of the graviton is equal to ^Iraviton^li^^^'graviton)^ ■ Therefore, in N=4 SYM all particles 

have masses, and the graviton is the lightest 
one. In such a theory, the large b dependence 
is determined by the mass of the lightest parti- 
cle dol, namely Aisi{s,b) exp {-mgraviton b) . 
This fact changes completely the ultra high en- 
ergy behavior, that has been considered in the 
previous section. Assuming that the graviton 
mass is small, we can distinguish four different 
kinematic ranges of energy in the case if Ra < 
< 1; Z^QnS > 1, but bn oc 


























1) 


\ 







Figure 7: The graviton(Pomeron) trajectory in N=4 SYM 
as it follows from Ref. iHl. 



z'^s 



< 



R\; and R\ 



< 



6q oc 



and 6q cx 



z'^s 



> 



, but bl 

1 / TTl ^ 

/ graviton 



s < 

■graviton- Nevertheless, we 
believe that the mass of the graviton should be 



^In the previous sections we called a'gravUon Just a'. 
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such that Ra ^ ^/mgraviton if N=4 SYM pretends to describe the main features of the strong inter- 
action. Indeed, we know experimentahy that the Ughtest hadron is the tt meson, and the large b de- 
pendence of the amphtude is proportional to exp {—b/lm-,^). For a massive graviton the amplitude falls as 
exp {—b/m graviton)- Therefore, to avoid contradiction with experiment, we need to assume that mgravUon > 
2m7r- Having this in mind, we will consider a modification to our formulae of the previous sections for 
the Glauber - Gribov approach in the case of the reggeized graviton, in three kinematic regions, which 
are z'^ gQ s < 1, where we can restrict ourselves to the exchange of one graviton in the nucleon scattering 
amplitude; z'^ g^ s > 1 but 6q oc z'^s < l/mgravUon < ( ™ ^^^^ region the multi graviton exchange 
could be essential); and the asymptotic region where z^ g^s > 1 but bQ (x z^ s > l/mgraviton < Ra- 
Of course, we can consider the kinematic region where {l/mgraviton Ins) > Ra-, but in this region nuclei 
behave in the same way as the nucleons, and we are not interested in this region. 



5.1 



< 1 



In this kinematic region we can restrict ourselves to one reggeized graviton (Pomeron) exchange, and 
use Eq. ( |5.5D instead of Eq. ( |2.32 ). Therefore, we have 



F2(g^x = QVs) = Ca'Q^ j cfb j dz z^ {Kf (Q z) + (Q z)) 2 Re {I 



(5.8) 



exp ( i glNc / dz \^{z')\'^ 2 z z ( cot ^ + ij F""^ J —^=-^ exp 



21ns 



Sib) 



Two features of Eq. ( |5.8p are quite different from Eq. (3.4), namely that the nucleon amplitude has an 
imaginary part and shows a different z dependence. Roughly speaking, in Eq. (^.§|), oc z"^ instead of 

(X z^ in Eq. (3^). The integral over z in Eq. (|5.8| ) can be evaluated using the steepest descent method. 
Using the asymptotic expression for the modified Bessel function, we reduce Eq. ( ^.SD to the following 
expression 



F2{Q^x = Q^/s) = (5.9) 
Ca'Q^ j dH j z^ dz e''^ ^ Re {1 - exp {ig^ N^^ijo) s^°'Hz z')^'' S (b) E {ln{z' .z)))} 

where E{ln{z'.z)) = W ^ ^ ^ exp | - ^ ^^^^^^ j and ^(j) = i cot ^ - 1 



Actually in Eq. ( |5.9| ), we need to integrate over z' , but we assume that the typical z' ~ 1/A, where A 
is a scale of hadrons (glueballs) in N=4 SYM, and we can replace it with some average value. 

It is convenient to introduce new dimensionless variables Q = Q z' , s = sz''^,S (6) = z''^S (6), z = z/z' , 
for which the equation for the saddle point reads as follows 
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Injl/zsp) 
Ins + In zsp 



E (ln(l/%p)) 



Rewriting Eq. (5.10) in the form 



In 



' - (jo - 1) i - ^ + w {zsp) 



where w is a smooth function of zsp, and t = In (s zsp)- The approximate solution for t is 



(5.10) 



(5.11) 



lns-^ln(^e-(0)l 
VX \S{b) J 

exp(^-^lnf^e-"'W 



-1+jo 

-1+io 



In s ; 



In s 



^ VX \S{b) J V ^ 

Using Eq. (|1|), we find that the DIS structure function is proportional to 



(5.12) 



F2{Q^x) oc expi-^-^Qzsp) = exp ( - (Q/Q,(x; ^))' V ^'^^ 



Q^exp 



V 



Jo 



(5.13) 



Q 



Q 



jo \glN,ajo)S{b) 



X s 



■ _ VX ln(l/zsp) \ V\ 
lns+ln%p 



In Eq. (|]13D we chose ^^^p, since it gives a larger contribution. The saturation momentum is equal to 



Qs {x\A) 



Jo 



Jo 



1 [x 



goNc^h) S (6) 



1 



_ ln(l/i5p) 
■JO In s-\-\n ZSP 



TT \/A{ln s+ln zgp) 



1 



1 



(5.14) 



where 



with j(A) = 1 



1 



jo - 1 



Vx \ Vx 

^From Eq. ( p. 14 ) one can see that F2 has a geometrical scaling behavior, if we neglect the log dependence 



of the saturation scale. The most interesting result is the fact that Qs oc (S'(6))^J(^' {l/x)^ ^ , 
At very large A, the saturation momentum is constant and does not depend on A and x. However, the 
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A dependance is more suppressed, namely A^^^'^^\ while the x dependence has a suppression, which is 
however a much weaker one (1/x)'^ ^ * . Such a behavior is similar to what we expect in high density QCD 
for the running QCD coupling |]3T| . 

5.2 Z^g^S > 1 but OC z"^ S < l/mgraviton < Ra 

In this kinematic region we have to take into account the multi-graviton interaction in the nucleon 
scattering amplitude. At high energy, the exchange of one Pomeron (reggeized graviton) induces an imag- 
inary part of the amplitude, as has been discussed, which increases with energy. Such an increase leads to 
a nucleon cross section of the order of 27r6o (x) , where 60 can be estimated using the following equation 

A%is,z,bo) « 1/2 (5.15) 



The nucleon amplitude for single reggeized graviton exchange can be evaluated using Eq. ( |5.2| ) and the 
fact that u{b) — > 5^/(2zz') at large b. Repeating the same procedure, we obtain that (with b = b/z') 

A^s^zM . ^(.)^° z^o - 2(lng + lni) ' ^'-^'^ 



iFrom Eq. (|5T|) and Eq. (|1§), we see that 6q ~ sz^. Therefore, 



F2 oc exp (^-gi - Const S (6) s^O"^ f^") (5.17) 
which leads to small values of the typical z = zsp, namely, 

(q/5(6))^ 

ZSP = (5.18) 

Jos 

At high energies, zsp is small, and the nucleon amplitude turns out to be small even at small values 
of b. Therefore, we do not need to discuss this region separately, and the answer is just the same as in the 
previous section. 

5.3 z'^g^S > 1 but 6^ oc S > l/nigraviton < Ra 



At such large impact parameters, we cannot use Eq. ( |5.lD and Eq. ( |5.5| ). The main contribution in 
this region stems from the exchange of the lightest hadron (in our case of the graviton) |^0| , which has the 
form given in Eq. ( 2.15| ), and can be written as 



A{s, b':^ z) oc ig^sz exp {-mgraviton b) (5.19) 



The typical impact parameter can be found from the equation A{Eq. ( 5.19 ); s, b) « 1/2, which gives 
60 = O-lfn graviton) In (So-^^s). Therefore for F2 we have 
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(Q^x = QVs) = (5-20) 

Ca'Q' I d% I dz e-'^' Re{l- exp {ig^ N, ^Jo) 2 vr S (6) (l/m^„,,,„J In^ {g^ )) } 

In Eq. ( pop , the main contribution stems from z (x 1/Q, which leads to 

F2{Q\x = Q^/s) = (5.21) 
Ca'Q' J d^ReS^l - exp (^i^^ iVc ^(jo) 2 vr 5 (6) (l/m^,,^,^ J In^ (^5o' ^ 



One can see, that F2 oc a' Q [In In \^gQ j j . However, such a behavior is vahd only in the 

restricted kinematic region when {l/m graviton) In {qq q^^) < ^A- For higher energies, we loose all the 
specifications related to the nucleus, and the nucleus interacts as a proton would do, but with the coupling 
constant g^ Nc A. 

6. DIS with nuclei: dipole model. 

In QCD, the DIS cross section can be written as a product of two factors |27, 15], namely the probability 
to find a dipole in the virtual photon, and the scattering amplitude of the dipole with the target. In this 
way the DIS cross section is given by the expression 

atot [DIS- Q\x) = I ^-^dH\^{Q;r,C)\'N{r,b,x) (6.1) 

where is the imaginary part of the scattering amplitude of the dipole with size r off the target, and C 
is the fraction of the energy carried by the quark of the dipole. We can try to generalize this equation to 
N=4 SYM, namely, 

d'b \^ {Q; r, z,C)\' Na {z, b, x) (6.2) 



The factorization of Eq. ( |0| ) is valid on very general grounds (see Ref. ^^), and should hold in any 
reasonable theory, since it is based on the structure of the interaction in time. In Eq. ( |6.2| ), we use the fact 
that the interaction due to the graviton exchange does not depend on the size of the interacting particles, 
(see Eq. (|2l5|)). We do not see any specific features for the dipole - target interaction, and thus we should 



be able to use for the nucleus amplitude {Na) the formulae that we have discussed in the previous sections. 
On the boundary, ^'(Q;r, C) is known and it is proportional to Ki [Qr), or to Kq (Qr), for different 
polarizations of the virtual photon with = Q'^C{1 — C)- We can reconstruct ^' {Q; r, C) using the Witten 



formula |32|, namely, 
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r(A) 

7rr(A - 1) 



j2 / 

a r 



^2 _|_ — -py 



J j Ko{Qr') with A± = -[d ± V(i2 + 4 



(6.3) 



Using the formulae 3.198, 6.532(4), 6.565(4) and 6.566(2) from the Gradstein and Ryzhik Tables, 
Ref. [33 1, and using the Feynman parameter (t), we can rewrite Eq. ( |6.3| ) in the form 



7rr(A- 1) 



(6.4) 

A+l 



r(A +1) 

7rAr(A- 1) 
1 

7r2A-ir(A- 1) ^ 7o + (1 - 

Using Eq. (|6.4D, we can rewrite Eq. (|6.1|) in the form 



I dt 



tz"^ + r'^t{l-t) + ^ 
X A-l 



atot {DIS;Q^ 



27r 



dz 



^2^-ir(A-l)^ io "^^^^ + (1-*)^'^ 



A-l 



(z2 + (l-t)r2)) 1^ i?e (^1 - expl^i^^^ ^z^S{b)j j (6.5) 

where we used the simple exchange of the graviton as in Eq. (^^). Using the asymptotical expression 
for the modified Bessel function, we can do the integral over z in saddle point approximation, and the 
equation for the saddle point zsp has the following form 



which leads to 



2 

E 

i=l 



\JtiiZgp 



=zsP + iglN,^ zip S {b) = 

2\ X 



(6.6) 



ZSP 



2 iQU 



\ glN^^S{h) \lQS{b)r 



<C r 



(6.7) 



In Eq. 
Prom Eq. ( 



^),Eq. ( |6^ ) and Eq. (|6.7|) , we introduced two variables, ti and t2, to describe |^ {Q;rX) 
and Eq. (|6.5D, we obtain 
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7rr2(A) fzspQ 



12\ 2A-2 ^3 
r2 y Q 



exp 



i=l 



1 , .5o'A^c 



Introducing a new variable Q = Q Yli=i VU, we can integrate over r using the steepest decent method. 
The equation for the saddle point reads as follows 



I X 



rsp 



2glN,S{b) Q Q 



(6.9) 



and 



the second term in Eq. 



rspd{r - rsp) dQ 2, 



2tt 



(zspQ^ 



2x 2A-2 , 



' S_P 

Q 



exp 



.2x1 



(6.10) 



with 



Ql 



2glN,S{b) 



oc 



yl3 1 



(6.11) 



Eq. ( |6.1Cl| ) shows geometrical scaling behavior, at least to within exponential accuracy. The saturation 



momentum of Eq. ( |6.11| ), has expected from the high density QCD A dependence, increases in the region 
of low X in the same way as for the DIS case, with the TZ current given by Eq. (|3.7|) . In general, Eq. ( |6.10D 
displays the same features as Eq. (|3.5D, (see also Eq. (|3.g|)). 



7. DIS in a shock wave approximation. 



The approach developed above, has to be compared with Ref.[26|, in which DIS with a nucleus target was 
considered in the framework of the shock wave approximation. In this paper the usual decomposition of 
the DIS cross section into two factors given by Eq. (^)[^, 28, |l^, which are the probability to find a 
dipole in the virtual photon, and the amplitude of the scattering of the dipole with the target, is used (see 
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Eq. ( |6.1D ) where C is the fraction of energy carried by the quark of the dipole. In Ref. |26] it is suggested 
to study the dipole-target amphtude in the semiclassical hmit of the dipole scattering, in the presence of 



the shock wave that was produced by the nucleus, in the spirit of Ref. ||3^. In this approach, the dipole is 
located at the boundary of the AdS^ space, and the two-dimensional surface of the string is characterized 
by X^^ = X^(r, a), (with /x = 0, . . . , 4), which depends on two coordinates (r, a). The string Nambu-Goto 
action takes the following form 



Sng = J dadrC = J drda-s/^detG^ where Ga,i3 = 9iJ.,v{X) daX^ dj^X'' , 



a, (3 



(7.1) 



In the presence of the heavy nucleus, the free metric of Eq. (2.18) has to be altered in order to take into 
account the energy-momentum tensor that describes the interaction of the dipole string with the nucleus. 
The modified metric is given by 



ds' 



— {^—2dxj^dx^ + {dx±)'^ + dz^^ + T 5{x-)dx-dx- 



(7.2) 



In Eq. ( |7.2[ ), we denote x± = '^"^'^ , where xq is the time in the normal four dimensional space. 



X4 = z. In Ref. |Q, T fj,z 6{x^ , suggested in Ref.[35|, is used. Using this assumption of Ref. |2C], the 

metric reduces to the expression 



ds' 



ds' 



ds' 



L 



— (2dx+dx- + IJ. z S{x^) dx_ + {dx±) + dz ^ 



z 

L2 



^ (2dx+dx^ + ^e(x_) e(a 



z2 V 
z2 



{ 



2a 



dt^ + 



1 + ^/ 
2a 



dx'i + {dxxf + dz' 



{dxsf + {dxi_f + dz^l 



(7.3) 



where a is chosen such that ^/2a = s^, and a ~ 2RaA/p+ cx (see Ref. ^ for details). In the 

last line of this equation, we omit the theta functions, since we are looking for the solution which does not 
depend on time (static solution [|26 



The static approximation is not well justified (see Ref.|39], which 
appeared after the first version of this paper we put on the net). However, the exchange of gravito, which 
interacts with the energy- momentum tensor (see Eq. ( |2.15D ) and which is responsible for the mediation 
of the gravitational force, is taken into account in this approximation. As we mentioned above, the main 
goal of this section is to confront the Glauber-Gribov approach for dipole -nucleus scattering, described 
in the previous section, with the static solution in the SW approximation. Although at first glance the 



solution of Ref. |26| does not reproduce the result of the Glauber-Gribov approach (see below), we will 
argue, that by changing the form of the Lagrangian of the string interaction with the nucleus, we are able 
to reproduce the Glauber-Gribov formula, in the static solution. Therefore, although it is plausible that 
one can learn some physics from the static solution, we believe however, that we can learn no more than 
is already derived from the Glauber-Gribov approach. 
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Using Eq. ( |7. 1| ) , Eq. ( [7. 3D and the following parameterization of X^, namely = t, X^ = x, X'^ = 0, 
X^ = and X'^ = z{x) as in ref. [26|, the action S is found to be equal to 



raV2 /•r/2 . a/A 1 

:,From the Euler-Lagrange equation, which has the form (for the static solution) 



S = dt f'^ dx C'''''^ with = — ^ Jil + z'^) (l-s2z4) (7.4) 



dx dz' dz 
as in ref. |26], the following solution is found 



a a r static a r static 

""^ = 0, (7.5) 



^/Aa fcgr^ 2 E^ (1/4) 1 

7rCoV2 U;^ax ^max J (27r)'^/^ 

while -Zrraoa; is the solution to the equation 



Cor = Zmax \/l " ^^ax (7-7) 



The amplitude in Eq. (|6^) is equal to |2§| 



N{r,x) = i2e{l -exp(i5(/u))} (7.8) 
Using Eq. ( |6.lD , the cross section for DIS has the form 

atot{DIS-Q\x) = I ^-^\^i,{Q;r,C)\'N{r,x) oc J {Qr) N {r,x) (7.9) 



In Eq. ( |7.9D , we omitted the integration over impact parameter, since in this simplified string approach 
we consider that the nucleus has the infinite extension in the transverse plane. As we have discussed 
above, such a simplified approach to the impact parameter dependence could cause a lot of difficulties, 
since DIS cross sections depend on the impact parameter distribution both in the nucleus and in the 
nucleon amplitude (see section 5). In Eq. ( [7.9| ) we simplified the wave function of the photon, which is 
known, by replacing it by Kq, since at large values of = <5^C(1 ~ 0) both components of the photon 
wave function for transverse and longitudinal polarized photons have the same behavior exp (— 2Qr). 

We expect that in DIS the typical r will be small, and therefore we try to find the solution to Eq. dTJ 



for which m = Cg <C 1. In this case, in ref.[^] three solutions have been found, which correspond to 
three different Riemann sheets of the cubic root, and which can be characterized by the index n = 0, 1,2. 
They are 
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m— »0 



1/ y/s for n = 0; 



i/ ^/s for n = 1; 



(7.10) 



Co r for n = 2; 



The solution with n = 2 is the only one that matches the Maldacena result |3€], for which m ^ 0. For 



this solution, we can take the integral over r in the second term of Eq. (7^) in Eq. ( [7.9| ) by the steepest 
decent method, with the saddle point 



rsP 



2^/2ttQ 



(7.11) 



One can see that 



msp 



4 4 2 



A \2 _2 \2 42/3 

" - ^ ^ « 1 for Q2 » x^A'l' 



oc 



and, therefore, in the kinematic region ^ A^^^, the second term of Eq. (|7.8| ) leads to an approach 
of the unitarity bound for the DIS cross section of the following form 



atot{DIS;Q^x) = J \^ {Q;r, C) \' - an (7.13) 

an = exp|-f%Vl (7.14) 



where the pre-exponential factor can be easily calculated. The saturation momentum Qg has the form 

2^/27r 



At large values of Q, Eq. ( 7.15| ) leads to a term of the order of xA^^^X/Qg, which corresponds to the 



twist expansion, with the anomalous dimension 7 = 1/2. The A dependence is in accordance with this as 
well [^], but the x dependence looks strange. — > at x — > 0, and therefore the theory predicts that for 
low X and > ^2/3^ ^j^g j^jg 

cross section is very small. 
It turns out that in the kinematic region <C A^A^/^, the n = solutions gives the largest 



contribution. Indeed, inserting this solution in Eq. (7^), one can find the saddle point in the integration 
over r, which is equal to 
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Evaluating m = Cq r^p s , namely 

one can see that in the region where <^ A^A^/^, we are deahng with small values of m, and we can 
use the solution of Eq. ( 7.1C| ). Then an in this case is proportional to 



an oc exp(-^-J with = 4 ^^^^^ _ oc ^ (7.18) 



The saturation momentum in Eq. ( [7.18 ) displays all the typical properties that we expect from high 
density QCD. 

It is worthwhile mentioning, that the solution with n = 1 leads to an oc exp i with the same 
saturation momentum, and it can be selected out since a should be positive. 

Both Eq. ( |7l^ ) and Eq. ( [tIsD have in common the fact that zf^ax^'^ turns out to be much smaller than 
unity ( -^max*^ It means that in the general equation for the action of Eq. ([t!^), we can consider 

z^s^ to be small, and we expand the action with respect to this parameter. In this case the contribution 
at high energy can be reduced to the following action 

geikonal ^ VXo^ f'''^ \/l " z'^ = CoUStsA^''^ f'^ dx Z^ \/ \ + z'2 (7.19) 

V27r J-rl2 J-r/2 

This action is closely related to the eikonal formula, as one can see from the second term of Eq. ( [7.19 ). 
Solving the Euler-Lagrange equation of Eq. ([7. 5]), we find that 



which leads to 



1 + = bnfL (7.20) 
z* 



zma. = i^^M%r, (^/2) (7.21) 
0Fr(3/4) 



Evaluating the integration over x in Eq. (7.19), we obtain the scattering amplitude in the form 



Nir,,) = ft {l - e.p (-C».,MV3 (-g;|^)- 

= i?e {l - exp (^-K^^/^sr^)} (7.22) 
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where we have absorbed all constant factors in the factor k. It is easy to see that Eq. (7.22) leads to 



an Of exp 




(7.23) 



with Qs given by Eq. ( 7.18 ). The difference between Eq. ( 7.23| ) and Eq. ( [7.18[) , as well as the fact that 
Eq. ( [7.14| ) does not hold, requires explanation. Referring back to Eq. ( |7.19| ), one can see that implicitly in 
geikonai ^ we neglected the part of the action of Eq. (|7.4| ) which does not depend on s. Since this contribution 
contains a factor of a oc 1/s in front, we can expect that this contribution is negligible at high energy. 
However, the integral over x can be divergent and compensates this smallness. In Ref. [26|, it was suggested 
that a subtraction in the action would cancel the divergence at 2 — > 0. The eikonal formula suggests a 
different type of remedy for this divergence, namely to introduce the action in the following way (compare 
with Eq. (O)) 



S = (7.24) 

dt / dxAr*"*^" = / dt dx AC'"''^ with AC'"''' = C{T^^) - C{T^u = Q) 

-00 J-r/2 Jo J-r/2 

which leads to 

^^static ^ J_^(l + ^/2)0(l_s2^4) _ A (7.25) 



Eq. ( 7.24 ) has a simple meaning, which is that we need to subtract the term which is responsible for 
the movement of the string in empty space during the period of time of the interaction, from the interaction 
induced by the energy-momentum tensor of the nucleus. One can see that the solution with the action 



given by Eq. ( [7.24 ) and Eq. ( 7.25| ), reproduces Eq. ( 7.23]) (see appendix) . 



Therefore, we can conclude that the shock wave approximation can be reproduced by the eikonal 
formula. It should be stressed that the eikonal formula is more general, since in the framework of this 
approach we are able to introduce the impact parameter dependence as well as the quantum corrections 
related to the reggeized graviton (Pomeron, see section 5). It is worthwhile mentioning that in Eq. ( [7. 181) 
the shock wave approximation leads to the same amplitude as Eq. ( |6.1[1| ), in the dipole approach. However, 
it should be stressed that the main result of Ref.|26|, that the dipole amplitude at high energy has a form 



N{r) oc l-exp-'"'3= with oc Const{x)A^/^ (7.26) 



holds in the approach with the action given by Eq. ( [7.24 ). Indeed, this result does not depend on the 



modification of the Lagrangian since for ^ 1 the action is the same in both approaches , namely, 

raV-2 _^ fr/2 ^ ^ 
-r/2 



l-aV2 /■r/2 /T 

S = i dt dx —s + ^ i Const{x) A^^^ r (7.27) 

Jo J-r/2 27r 
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The last equation comes from the equation of motion which leads to z' = 0. This saturation momentum 
Qsi-A) oc Const{x)A^^^ needs an explanation since it does not appear in the Glauber-Gribov approach. 
First, the contribution of Eq. ( 7.26| ) to the DIS cross section (see Eq. (|6.l|)) has the form 

For Q > Qs{A.) a [DIS;Q^,x) vr/Q^ which corresponds to 1 in Eq. ( 7.26 ). If we replace Kq by its 
asymptotic behavior a oc (1/Q^) x (1/(Q + Qs{A)). The physical meaning of Qs{A) is rather obvious: 
during the passage of the dipole through the nucleus the transverse momentum (Q) can get an additional 
momentum AQ due to elastic rescattering with the nucleons , namely 

AQ (X X number of colhsions = ^^^^ (7.29) 

Rn 

where = l/Rjy is the typical transverse momentum for elastic scattering with one nucleon {R^ is the 
nucleon radius). In the Glauber-Gribov approach, however, q^ oc 1/Ra due to the nucleus form 
factor (see Eq. ( |2.8| )). In the shock wave approach the nucleus wave function has not been taken into 
account and nucleons can have unrestricted transverse momenta. Therefore, we consider this momentum 
as the artifact of the shock wave model in which, we believe, we need to specify the DIS as scattering with 
Q > AQ ~ Qs{A) if we are interested in finding the total cross section. However, this Qs{-A) can manifest 
itself in the inclusive production leading to the situation with two characteristic momenta that has been 
advocated in Ref. ||3^. It should be stressed that Eq. ( 7.2S| ) is written for a string with the fixed transverse 
coordinate (see Eq. ( [7. 4]) = x, X'^ = 0). In the general case AQ^ = -^A^^^. The second comment 
on Eq. ( 7.29D is that we considered the rescatterings which are instantaneous in accordance with the static 
solution. In the region Q > Qs{A), the contribution, given by Eq. (7.28), is small and the value of the 
total cross section for DIS is determined by the saddle point approximation (see Eq. ( f7.23D ) which is the 
same both in the shock wave approximation and in the Glauber-Gribov approach. 



8. Conclusions 



It is our common wisdom nowadays that N=4 SYM , which can be solved at large coupling values, can 
provide us with some knowledge of what potentially lies in the confinement region of QCD. However, the 
first analysis of high energy DIS scattering, performed in Refs. [^, |lO|, demonstrated that the high energy 
scattering in N=4 SYM looks quite different from what has been known so far. Contrary to the usual 
expectations based on perturbative QCD and the parton model, that the main process at high energies is 
multiparticle production, it was found in Refs. g [H that in N=4 SYM the major contribution originates 
from quasi-elastic scattering. This also contradicts what is known from data. 

The goal of this paper was to develop the Glauber-Gribov description of DIS on a nuclear target within 
the N=4 SYM, which should help to see the key features of high energy scattering in a more transparent 
way. For this purpose we employed the eikonal approximation which has been developed for N=4 SYM in 
Refs-I^, ^, ^, ^, H^. Our results can be summarized as follows. 
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1. We derived the Glauber-Gribov formula (see Eq. ( p. 41 ) and Eq. (^?5|)), and showed that for the case of 
graviton exchange, this formula displays the same general properties, such as the geometrical scaling 
behavior, as in the case of the high density QCD approach. 



2. We demonstrated that graviton exchange indeed leads to a total cross section which is dominated 
by quasi-elastic re-scatterings. However, we found that the quantum effects responsible for graviton 
reggeization give rise to an imaginary part of the nucleon amplitude. This imaginary part, enhanced 
by multiple interactions, results in a DIS which looks similar to one predicted by the high density 
QCD, (see Fig. |). 



3. We concluded that in N=4 SYM the impact parameter dependence of the amplitude is essential, and 
the entire kinematic region can be divided into three regions. In the first region [z^ g^Nc < 1), we 
can use the eikonal formula with a single graviton or reggeized graviton exchange for the nucleon 
amplitude. In the second kinematic region, g^Nc > 1 but 69 ^'^^ < ^1 "'^graviton < -^A' ^'^^ 
multi-graviton exchange in the nucleon amplitude may become important. However, we found that 
this is not the case and still the single graviton exchange dominates. In the third kinematic region 
{z^ g^Nc > 1 and ^ / ni^graviton < (X z^s < R^), the multi-graviton exchanges in the nucleon 
amplitude must be included, and the related modification to the amplitude are discussed in section 
5.3. 



In this paper, we considered mostly the DIS of the TZ current with the target. However, in the last 
two sections, we discussed the traditional approach to DIS based on the factorization given by Eq. (|6.l|) . 
We considered DIS in two different ways. In the first one we generalized the usual dipole formula to N=4 
SYM. We derived the probability to find a dipole in the virtual photon, in AdS^ space, and considered for 
the dipole scattering amplitude the eikonal formula. In the second approach, we revisited the shock wave 



approximation that has been developed for DIS in Ref. |26], and we showed that in this formalism we can 



also use the Glauber-Gribov approach for DIS in the region of r ~ Q/A^^^^/s . However, the Glauber- 



Gribov approach suggests a different way to renormalize the interaction Lagrangian proposed in Ref.|2f:]. 



After such modification of the original formalism of Ref. |26|, both approaches, namely, the dipole model 



and the shock wave approximation give the same result for r ~ Q/{A^/^^/s). We gave the interpretation 
of the appearance of the new saturation momentum Qs{A) that does not depend on energy [^] and argue 
that in the shock wave approximation we should consider only DIS with Q > Qs{A). For such large 
values of Q the shock wave approximation with our modified Lagrangian reproduces the same result as the 
Glauber-Gribov approach. 

In general, we conclude that N=4 SYM does not lead to any obvious contradiction, either with the 
high density QCD, or with experimental data. Therefore, we hope to learn something valuable about the 
confinement region from the exact solution in N=4 SYM, relying on the AdS/CFT correspondence. 
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A. Shock wave approximation for DIS with our hypohesis on renormahzed Lagrangian 



In this appendix we consider the shock wave approximation to DIS suggested in Ref. |2£], with our 
hypothesis on the renormalised Lagrangian. As has been mentioned, we assume that the static AdS^ 
renormalised lagrangian is the regular AdS^ lagrangian with a nucleus present, minus the vacuum AdS^ 
lagrangian, where the nucleus is not present. The expression to such a renormalised lagrangian is given by 
the following expression 



^^6" - C (T^,) - C {T,, = 0) = - (A-1) 



V2X 1 



where = ^ {I + z'^) {1 - s^z^) and = ^ ^x/(l + z'2) (A-2) 

2tt ztt z'^ 



The Euler - Lagrange equation for C}^"^ takes the form; 



d z dx 



dz dz dx V d z' J dx \ dz 



+ =0 (A-3) 



The various terms appearing in Eq. ( [A- 3D can be calculated from Eq. (A-1), namely 



dC 



nuc 



dz 

and — (^^^\ 
dx \ dz' J 
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g£vac 
dz 
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_f £vac 

z 



z 



■C 



nuc 



(?x \^ 1 + 



l + z 



12 



{l + z' 



2^2 



./2 



^nuc ^ 



+ 



z'z" 



z''^z" 



./2 



similarly ^ J 



^vac ^ 



I + dz 

/2 ^£vac 



I + Z'2 (1 + ^/2)2^- ' l + dz 

Plugging Eq. ( [A^ ), Eq. ( [A^ ) and Eq. ( |A76t ) into Eq. ( [A^ ) gives the result 



(A-4) 
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S2z4 











Recall that one can express z" as (1/2) dz''^ /dz, hence Eq. (|A-7| ) simplifies to 



(A-7) 



2 l + z 
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1 (9z 



12 



2 5z 



Vl — S^Z 



2^4 



dz 



/2 



l + z 



12 



dz 



Z y/l- S^Z 



2-4 



(A-8) 



Integrating between z' (z) and z' (z^) = 0, where extremum point of the string, one arrives 

at the result 



J2 



, Zm . 



l + ^/l 



„2^4 



i + Vi 



s'^z^ 



(A-9) 
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^From Eq. ( [A-9| ), one can find that 



f 

Jo 



^/i (x - r/2) 



(A-10) 



where ^ = \/sz. In Eq. ( [A-10| ), the half of the string where z' > is chosen, and we integrated over 
X from —r/2 to x. The maximal value of ^ = can be found from the following equation 



/sr/2 



(A-11) 



We have not yet found the expression for the function H {^m, £,m) through known functions, but the 
figure of Fig. ^ and Fig. ^ demonstrates the behavior of this function. The key difference with the solution 
proposed in Ref. is the fact that H {^m,^m) of Eq. (|A-11 ) has only one solution in the region of small 



^fn, while H{^rn,^m) of Ref. |26] has two solutions (see Fig. g). We can simplify the integrand by its 
expression at low ^, namely, 



H 



Low ^ 



JO 



dCU 





Changing the integration variable to = sinO, then Eq. ( A-12| ) becomes; 



(A-12) 



H 



Low ^ 



2 ^ ./o VsinI 



(A-13) 



Finally changing the variable of integration once again to V sin 9 = t, Eq. ( A-13| ) reduces to 



H 



Low 



i/im dt 





s{x — r/2) 



ellpt{arcsin(^/^m) ,0} 



(A-14) 



where ellpt {(p, k) is the elliptic function defined as 



ellpt {4), k) 



dt 



Vl - kH'^Vl^ 

At large values of ^, expanding the integrand at large values of ^, we obtain 



(A-15) 
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2 + 2^1^ 



Sm S 



2 + 2^1^ 



^/s {x - r/2) (A-16) 



Fig. P and Fig. |9| show how the simpUfied equations (Eq. (|A-12| ) and Eq. ( A-16|) ), describe the exact 
function H (£„, given by Eq. (ITT^) . 
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Figure 8: Function H{^rn,£,m) versus Cm for small 
values of The solid line shows the function 

H{£,ra,Cm) given by Eq. ( A-1C1| ), and the dotted hne 
is the same function for the solution given in Ref. 
p6[ while dashed line describes the approximation 



of Eq. ( |A-12| ). 




Figure 9: Function (£m,£m) versus for large 
values of £m. The solid line shows the function 



) given by Eq. ( [A-IOD , dotted hne is the 
approximation by Eq. (A-16). 



Using Eq. ( |A-12| ), one can easily see that Eq. (|7.4| ) with Lagrangian of Eq. ( A-1 ) reproduces Eq. ( 7.23 ), 
which we obtain from the eikonal formula. 
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